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We show that the phase boundary between the paramagnetic metal and the nonmagnetic Mott in- 
sulator for the Hubbard model on a triangular lattice obtained by Yoshioka et al. in arXiv:0811.1575 
does not correctly represent that of the thermodynamic limit but is an artifact of the 6 by 6 lattice 
they rely on. After the system size extrapolation, the phase boundary is located at U/t ~ 5.2 as 
proposed by Morita et al, J. Phys. Soc. Jpn. 71 (2008) 2109 and in contrast to Yoshioka et al. 
Here, U is the onsite Coulomb repulsion and t is the nearest-neighbor transfer. 



Recently, Yoshioka, Koga and Kawakami up- 
loaded a paper entitled "Quantum Phase Transitions 
in the Hubbard Model on Triangular Lattice" on 
arXiv:0811.1575,[H They have discussed phase diagram 
of the Hubbard model on a triangular lattice in the 
ground state as a function of the onsite Coulomb re- 
pulsion U and the nearest-neighbor hopping t by using 
the path-integral renormalization group (PIRG) method. 
In particular, they have reproduced the existence of the 
nonmagnetic insulating phase near the Mott transition 
in agreement with the previous PIRG result by Morita 
et al. However, they found a phase boundary between 
a paramagnetic metal and a nonmagnetic insulator at 
U/t ~ 7.7 with a large jump of the double occupancy 
D = {ni^rii ^), AD ~ 0.07, quantitatively contradict- 
ing the r esult of Refi, where the phase boundary has 
been found at U/t ~ 5.2 with a small jump AD (less 
than 0.025) of the double occupancy suggesting a weak 
first-order or continuous Mott transition. Here rii^ is 
the number operator at site i with spin a. Yoshioka et 
al. have claimed that their "new" algorithm of itera- 
tive method allows going beyond the possible trap in the 
metastable state. According to them, this makes the dif- 
ference from the result by Morita et al. 

Here, we show that the difference does not result from 
the issue of the trappingbut from size effects, and the 
result by Yoshioka et aZ.[lJ is an artifact of their calcu- 
lation performed only for 6 x 6 lattice. We have indeed 
rechecked the data obtained in the calculation by Morita 
et al. and confirmed that the data for the 6x6 lattice are 
essentially consistent with the properties of the ground 
state by Yoshioka et al.. Yoshioka et al. claimed that 
their data have better reliability because they have im- 
plemented an improved procedure in the iteration of the 
path integral operation, where more independent paths 
of the Stratonovich- Hubbard variables are considered. In 
fact, although it has not explicitly been described, we 
have already implemented in the calculations presented 
in Ref. essentially the same algorithm of their itera- 
tive algorithm as a derivative one of our standard proce- 
dure of PIRG. This minor algorithmic extension was also 



supplemented by further implementation of the adaptive 
interval of the imaginary time slice etc. 

The ground state energy of the 6x6 lattice at U — 6.5 
and t — 1 is —0.791 ± 0.007 per site and the double 
occupancy D = (rt,;f rt,;j) is 0.153 ± 0.004 in our calcula- 
tion. In Fig. [TJ we show our data for momentum dis- 
tribution ?i(q). These are obtained after the quantum 
number projection of the total spin. [3] The results show 
good consistency with the results by Yoshioka et al. In 
addition, the metastable state shown by Yoshioka et al. 
in the range of 6 < U < 7 is also consietent with the 
metastable states we obtained, which have not been dis- 
cussed in the paper by Morita et al., because these are 
simply excited states. The data are also consistent at 
U larger than 8, where the nonmagnetic insulating state 
seems to be well stabilized. In the interval 7 < U < 8, 
our results show that only the state with a nonmagnetic 
and insulating character with relatively small D around 
0.1 is stabilized in contrast to the result by Yoshioka et 
al. In Fig'121 we show the U dependences of the energy E, 
double occupancy D and the spin structure factor 5(q) 
at the peak, both for the ground state and the excited 
state, showing the behaviors similar to those by Yosh- 
ioka et al. aside from the small discrepancy mentioned 
above. Here S'(q) is the Fourier transform of the spin 
correlation — «i,i)(?ij,T — n j,l))/^- The results for 

the 6x6 lattice have, therefore, overall consistency be- 
tween these two independent calculations except for a 
small discrepancy in the region 7 < U < 8. 

A problem of the small size system such as the 6x6 
lattice is that it is not clear what type of phase an ob- 
tained state belongs to. Actually, the ground state we 
obtained shows the momentum distribution illustrated 
in FigQ] which could be interpreted both by a metal and 
an insulator in the range of the interaction 6 < U < 8. In 
any case, for the 6x6 lattice, D and S(q) seem to show 
a sharp crossover or level crossing between U — 1 and 8. 
This transition was identified as the metal-insulator tran- 
sition by Yoshioka et al. However, as noted above, the 
distinction between metals and insulators is not definitely 
clear when one considers only one size of lattice, and the 
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FIG. 1: Momentum distribution n(q) of the ground state at 
U/t=6.5 for 6x6 triangular lattice at half filling. 




only the results up to the 6x6 lattice, it would not be 
possible to determine the phase boundaries. 

In addition, small-sized systems such as 4 x 4 and 6x6 
often show behaviors very different from larger size sys- 
tems. This was already cautioned by Kashima et aL[3j, 
and Mizusaki et aZ.[3], where a clear nonmagnetic insu- 
lating phase does not show up for the system size equal 
to or smaller than 6 x 6 for the square lattice Hubbard 
model with the next-nearest neighbor transfer, while it 
exists in the system larger than 6x6 and the size extrap- 
olation supports that it survives in the thermodynamic 
limit. This is equally true for the triangular lattice. In 
FigOl we see a sharp crossover (or level crossing) be- 
tween U — 7 and 8 for the 6x6 lattice. However, when 
the system size is increased, the sharp crossover for the 
6x6 lattice diappears already at 8 x 8 and the transi- 
tion is identified at much smaller U. In fact, the phase 
boundary becomes stable among these larger-sized sys- 
tems. The system size dependence of the charge gap A c 
defined by (E n+ 2+E n - 2 — 2E n )/2 shown in Fig. [3]clearly 
shows that the region 5.2 < U has a nonzero charge gap 
at half filling. Here E n is the ground state enegy of the n 
particle systems. Accordingly, the momentum distribu- 
tion shows insulator-type smooth feature without a jump 
at the Fermi leveljH as taken from the data in the same 
calculation with Morita et ai.pj shown in FigQJ It is 
unlikely that the obtained states for the larger-sized sys- 
tems are trapped in the metastable state, because essen- 
tially the same state is obtained after the expansion of 
the Hilbert space even when we start from different initial 
configurations whichever they are metallic or insulating 
Hartree-Fock states. During the process of the increasing 
basis states with iterations, the states converge to similar 
states. 



FIG. 2: U dependence of (a) total energy per site E (b) dou- 
ble occupancy D and (c) peak value of S'(q) for the ground 
state (red open circles ) and a typical excited state (blue open 
squares) of 6 x 6 lattice of triangular lattice at half filling. 




size scaling is imperative. Actually in our experience, it 
is hard to distinguish metals from insulators and antifcr- 
romagnetic states from paramagnetic ones from the data 
of a single system size. For example, if we see Figs. 8, 10, 
II and 13 in Refi for the square lattice Hubbard model 
with the next-nearest neighbor transfer, and if we knew 



FIG. 3: U dependence of charge gap in the thermodynamic 
limit (red filled circles). Open blue square is located at U = 
5.2 where the double occupancy shows a sharp change @] (see 
Fig.©. 

System size dependences of the double occupancy D for 
several different choices of U plotted in Fig. [5] reveal Q 
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U/t = 5.6 
t'/t = 1.0 




FIG. 4: Momentum distribution on 14 x 14 lattice at half- 
filling for Hubbard model on triangular lattice at U/t = 5.6 



of the charge gap A c , the momentum distribution n(q), 
and the spin structure factor S (q) , our conclusion for the 
phase boundary between the paramagnetic metal and the 
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that, though D stays above 0.16 up to U/t = 6.0 for the 
6x6 lattice, the extrapolated values in the thermody- 
namic limit continuously decrease in the insulating states 
above U/t > 5.2 even to D ~ 0.12 at U/t = 6, while D 
below U/t — 5.2 changes its system size dependence to 
an increasing function with increasing system sizes. This 
leads to D higher than 0.2 even just below the Mott tran- 
sition point U/t = 5.2 as we see in FigJSJ [2, 5] This high 
value of D in the metallic states near the Mott transi- 
tion has been already discussed by Tahara et al. in the 
case of the square lattice with the next nearest neigh- 
bor transfer. [6j The jump at the Mott transition appears 
to be small as we see in FigJBl suggesting a weakly first 
order or continuous transition. 
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FIG. 5: System size dependence of double occupancy D for 
the triangular lattice at half filling for several choices of U. 
Here the number of sites is N. 

From the system size dependence and U ft dependence 



FIG. 6: Double occupancy D in the thermodynamic limit 
after finite-size scaling for triangular lattice at half filling. 



nonmagnetic insulator identified around U = 5.2 with 
a weakly first-order character or continuous transition 
proposed by Morita et al. 0] is solid and we conclude that 
the results by Yoshioka et al do not represent the phase 
boundary of the thermodynamic limit. 

At larger U, the triangular lattice system may undergo 
a further transition to some type of antifcrromagnctic or- 
der. This possibility has not been seriously considered by 
Morita et aL0 and calls for further studies. Yoshioka et 
al. claimed the transition to 120°-type antiferromagnetic 
order at U ~ 10. However, this possible phase bound- 
ary has also to be examined with a careful analysis of 
finite size effects together with other possibility of an- 
tiferromagnetic order with longer-ranged period such as 
an up- up- down- down structure as proposed in the case of 
the square lattice Hubbard model with the next-nearest 
neighbor transfer. [4| This is an intriguing open question 
further to be studied. 
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